From chiral quark models to high-energy processes*** 
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We show the results of low-energy chiral quark models for soft matrix 
elements involving pions and photons. Such soft elements, upon convo- 
lution with the hard matrix elements, are relevant in various high-energy 
processes. We focus on quantities related to the generalized parton distri- 
butions of the pion: the parton distribution functions, the parton distri- 
bution amplitudes, and the generalized form factors. Wherever possible, 
the model predictions are confronted with the data or lattice simulations, 
where surprisingly good agreement is achieved. The QCD evolution from 
the low quark model scale up to the scale of the data is crucial for this 
agreement. 

PACS numbers: 12.38.Lg, 11.30, 12.38.-t 



Low-energy chiral quark models are designed to describe the the dynam- 
ics of the pion and its interactions. Their treatment at the large-iV c level 
is particularly straightforward, as it requires a rather simple task of com- 
puting the one-quark-loop diagrams with insertions of pion or gauge-boson 
vertices. Despite numerous works [1-24] on the applications of low energy 
chiral quark models to high-energy processes, their usefulness and predic- 
tive power is not well known to the high-energy community. It has been 
widely believed that because the non-perturbative structure of the pion is 
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such a tremendously difficult problem, it is better to parametrize the matrix 
elements (for instance, the parton distribution function) by some suitable 
form. Then this form is inserted as the initial condition to the QCD evolu- 
tion equations, which bring it to the relevant and experimentally accessible 
scales. Clearly, that way one can only relate the quantities at different scales 
and probe the evolution itself, but no genuine non-perturbative dynamics 
present in the soft matrix elements is investigated. In the case of the pion 
its low energy properties are expected to be dominated by the spontaneous 
breakdown of the chiral symmetry, a feature which is implemented in chiral 
quark models but is largely ignored in many high energy studies. 

In our approach two basic elements are crucial: the low-energy dynami- 
cal quark model itself, as well as the QCD em evolution, moving the model 
predictions from the low quark model scale to higher scales, where the ex- 
perimental or lattice data are available. Without an operational definition 
of the low energy scale Qo and the QCD evolution relating Qo to Q one 
cannot compare the model predictions to the data at the scale Q. This talk 
is mainly based on results published recently in Refs. [25,26]. 

The theoretical framework is set by the Generalized Parton Distributions 
(GPDs) [27-35]. For the case of the pion, the GPD for the non-singlet 
channel is defined as (we omit the gluon gauge link operators, absent in the 
light-cone gauge) 



ezab H«> I=1 (x, (,t)=J -^e ix P +z -(Ap')\mi + m ^3 k») 

while in the singlet case 

< 5 ab ^ / = (x,C,t)=/^e^ +2 V(p , )l^(0)7 + ^(^)k a (p)) 



z+=0,z x =0 
(1) 



=0,z x =0 



(2) 

where the kinematics is set by p' = p + q, p 2 = m 2 , q 2 = —2p ■ q = t, and 
£ = q + /p + , which denotes the momentum transfer along the light cone. 
GPDs provide very rich information of the structure of hadrons, which may 
eventually come from such exclusive processes as ep — > ep^f, ^p — > pl + l~ , 
ep — > epl + l~ , or from the lattice calculations. 

Formal properties of GPDs are most elegantly written in the symmetric 
notation: £ = ^r^, X = fr0§> with < £ < 1, — 1 < X < 1, where one 
finds 

H I= \X, S, t) = -H I=0 (-X, £, t), H I= \X, S, t) = H^i-X, e, t). (3) 
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For X > one has 

H I=o ' 1 (X,0,0)=q s ' NS (X), (4) 

with q(x) 1 being the standard parton distribution functions (PDFs). The 
following sum rules hold: 

V£ : JdX H I=1 (X, £, t) = 2F v (t), 

JdXXH I=0fv e ^ - na - ( * ot2 ' 



\X,£,t) = 29 2 (t)-2ee l (t), 



where Fy(t) is the electromagnetic form factor, while 9\{t) and 9 2 (t) are 
the gravitational form factors [36], related, correspondingly, to the charge 
conservation and to the momentum sum rule in the deep inelastic scattering. 

More generally, the polynomiality conditions [27], following from the 
Lorentz invariance, time reversal, and hermiticity, state that 



/ alX X^ H I=1 (X, C, t) = 2 ]T A 2Wl (t)e\ 
JdXX^H^X^t) = 2£A 2j+2)2i (t)£ 2i , 



where A's are the generalized form factors ( GFFs). These form factors can 
be also expressed as matrix elements of the form 

(ir + (j/)\ «(0) 7 {M iD^iDf* . . . W^-^ u(0) \tt + ( P )) = 
2jP {MpMi . . . P^-i} A n0 (t) + 

n 

2 •••<f k - 1 P tik ...P^ n - l] 2- k A nk {t), (5) 

fe=2,4,... 

where P = p' + p, hence the GPDs may be viewed as an infinite collection 
of GFFs. The positivity bound [37, 38] states that 



\H q (X,£,t)\ < ^q(x iQ )q(x out ), £ < X < 1, (6) 

where x- m = (x + £)/(l + £), x ou t = (x — £)/(l — 0- Finahy, a low-energy 
theorem [39] 

tf /=1 (2z- 1,1,0) = 0(z) (7) 
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Fig. 1. Left: model prediction for the valence parton distribution in the pion evolved 
to the scale of Q = 4 GeV (band). The width of the band indicates the uncertainty 
in the initial scale Qo, Eq. ((9|. The data points come from the Drell-Yan E615 
experiment [40]. Right: same, evolved to the scale Q — 0.35 GeV. The data come 
from the transverse lattice calculations of Ref. [41] 

holds, relating the GPD to the pion distribution amplitude (DA), denoted 
as 4>. 

We stress that the above-listed relations and bounds impose severe con- 
straints on the possible form of the GPDs. All these formal requirements 
are naturally satisfied in our quark-model calculations [25]. 

With C = t = 0, the GPD becomes the usual PDF, cf. Eq. (|4|). Long 
ago Davidson and one of us (ERA) found that in the Nambu-Jona-Lasinio 
model [1] q{x) = 1. This result pertains to the low-energy chiral-quark- 
model scale, which is a priori not known and has to be determined. At 
this scale the gluons have been integrated out and effective quark degrees of 
freedom are the only degrees of freedom. Thus, all observables are saturated 
by the quark contribution, in particular, the momentum sum rule. From 
experiment, the momentum fraction carried by the valence quarks is [42,43] 

(x) v = 0.47(2) at Q 2 = 4 GeV 2 . (8) 

We evolve this value backward with the LO DGLAP equations down to 
the scale where quarks carry all the momentum, (x) v = 1. This procedure, 
assuming the matching of the quark model results and QCD at the quark 
model scale, yields the quark model scale 

Qo = 3L3±fg MeV, (9) 

where the range reflects the uncertainty in 0. 
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Fig. 2. Left: model prediction for the pion distribution amplitude evolved to the 
scale Q = 2 GeV (band) and compared to the E791 di-jet measurement [48]. The 
width of the band indicates the uncertainty in the initial scale Qq, Eq. ©. The 
dashed line shows the asymptotic form 4>{x, oo) = 6x(l — x). Right: same evolved 
to the scale Q — 0.5GeV compared to the transverse lattice data [41]. 



We admit that such a low scale, where ev(Qo)/(27r) = 0.34, makes the 
evolution very fast for the scales close to the initial value Qq. This calls for 
improvement of the evolution. However, the NLO effects are small [9]. On 
the other hand, the removal of the Landau pole by means of the analytic 
QCD method [44-47] results in a too slow evolution, incapable of reconciling 
the quark-model results with experimental data. Thus in our studies we 
employ the simplest LO DGLAP QCD evolution, baring in mind that it is 
good at intermediate values of x. 

The results of our procedure for the non-singlet PDF (valence quark 
distribution) of the pion are shown on the left side of Fig. [TJ We have evolved 
the quark model result ([5]) from the scale Q up to the scale Q = 4 GeV 
corresponding to the E615 experiment [40]. We note a very reasonable 
agreement for all the range where the data are available. On the right side 
we show the same quantity evolved to the rather low scale of Q = 350 MeV 
and confronted with the transverse lattice data [41]. Again, the agreement 
is quite remarkable, showing the features of transverse-lattice calculation, 
designed to work at low-energy scales. 

Next, we look at the DA, related to the GPD via the sum rule ([7]). Here 
the evolution is carried out with the LO ERBL equations. The results are 
displayed in Fig. [21 again in fair agreement with the data, especially for the 
lattice case on the right-hand side. Thus the GPDs for the special kinematic 
cases of Eqs. (|4l7p are well reproduced in chiral quark models supplied with 
evolution. 

In Ref. [25] we provide expressions for the full form of the GPDs in the 
NJL model and in the Spectral Quark Model [50]. These expressions have 
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Fig. 3. The electromagnetic form factor (left) and the quark part of the gravita- 
tional form factor, A 2 o (right) in the Spectral Quark Model compared to the lattice 
data from Ref . [49] . The band around the model curve indicates the uncertainty in 
the quark momentum fraction and the m p parameter in the Spectral Quark Model. 



a non-trivial structure, in particular they do not exhibit factorization in the 
t and x variables, while satisfying all formal requirements listed above, in- 
cluding polynomiality. Since there is no data for the full kinematic range 
for the GPDs, we display the results for the generalized form factors. Inter- 
estingly, there is recent information on these objects from the full QCD on 
the lattice [49,51]. The Spectral Quark Model vector form factor and the 
quark part of the gravitational form factor of the pion, A20, are compared 
to these lattice data in Fig. [3J We note a very good agreement. In the 
Spectral Quark Model the expressions are particularly simple, 

F S y QM (t) = -£L (10) 

m~ — t 

d SQM {t) = e S Q M {t) = H£ log ( . 

t \m 2 p -t) 

We note the longer tail of the gravitational form factor in the momentum 
space, meaning a more compact distribution in the coordinate space. Ex- 
plicitly, in the considered chiral quark models we find the simple relation 

2(r 2 ) e = (r 2 )v (11) 

In Fig. U] we show our predictions for the higher-order generalized form 
factors, evolved to the scale of 2 GeV. In Ref. [52] we have pointed out that 
these moments evolve in a very simple way with the LO DGLAP-ERBL 
evolution. One finds triangular structures, for instance in the non-singlet 
odd-n case 



A w (t,Q) = LiA 10 (t,Qa), 
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A 32 (t, Q) = -{L x - L 3 )Aio(t, Qq) + L 3 A 32 (t, Q ), 
5 

A u (t, Q) = ^(VLi ~ UL 3 + 5L 5 )A 10 (t, Qq) 
+-(L 3 - L 5 )A 32 (t, Qq) + L 5 A 54 (t, Qq), 



A 30 (t,Q) = L 3 A 30 (t,Q ), 

A 52 (t, Q) = -(L 3 - L 5 )A 30 (t, Qq) + L 5 A 52 (t, Q ), 



A 50 (t,Q) = L 5 A 50 (t,Q 



(12) 



where 



Ln -{^m) ' (13) 

and 7i denotes anomalous dimensions, in the present case for the vector 
operator. Certainly, L\ = 1, and the vector form factor, corresponding to a 
conserved current, does not evolve. Similarly, the gravitational form factors 
do not evolve. The higher-order form factors, however, do change due to 
evolution. Similar structures of equations to (fT2l) appear for the even n, as 
well as for the singlet case [52]. At the moment, there is no data to compare 
to for the higher-order form factors, but our predictions of Fig. H] can be 
tested in future lattice studies. 

However, already at present there exists useful information at t = 0. 
Below we compare our model values of the higher-order form factors at 
t = to the lattice data provided in Sec. 7 of Ref. [49] . With the short-hand 
notation (x n ) = A n+ ifi(0) one finds that at the lattice scale of Q = 2 GeV 

(x) = 0.271 ±0.016, (14) 

(x 2 ) = 0.128 ±0.018, 

(x 3 ) = 0.074 ± 0.027. (lattice) 

while in the chiral quark models we get, after the LO DGLAP evolution to 
the lattice scale, 

(x) = 0.28 ±0.02, (15) 
(x 2 ) = 0.10 ±0.02, 

(x 3 ) = 0.06 ± 0.01. (chiral quark models) 
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Fig. 4. (Predictions for the higher-order generalized form factors ^3,24 and Aj^j of 
the pion from the Spectral Quark Model evolved to the scale Q = 2 GeV. 

The model error bars come from the uncertainty of the scale Qo in Eq. ([9]). 
As we can see, the agreement falls within the error bars. 

To briefly summarize the results presented in this talk, we state that 
the chiral quark models supplied with the QCD evolution work well for a 
wide variety of quantities related to the GPDs of the pion. Thus they may 
provide valuable insight into the non-perturbative dynamics behind the soft 
matrix elements. Further details can be found in Refs. [25,26]. 



Since at this meeting we are remembering Professor Jan Kwiecihski, I 
wish to end my contribution by showing one of the last photos of our dear 
teacher and friend sitting happily among his colleagues. The picture was 
taken on the occasion of handling to Jan the issue of Acta Physica Polonica 
B [53] (placed at the corner of the table) dedicated to him in honor of his 
65th birthday. The event took part in the Institute of Nuclear Physics in 
June 2003, just two and a half months before his premature passing away. 
We all miss you, Jan! 
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Fig. 5. 12 June 2003, left to right: Mariusz Sadzikowski, Wojciech Florkowski, 
Michal Praszalowicz, Kacper Zalewski, Krzysztof Golec-Biernat, Andrzej Bialas, 
Jan Kwiecinski, Maria and Wieslaw Czyz. 



